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We perform a nonperturbative tuning of the coefficients in the 0(a 2 ) improved action for staggered fermions. 
The mass splitting for the pions of different doubler flavor is used as a measure of the symmetry breaking effects 
introduced by 0(a 2 ) discretization errors. We find that the flavor nondegeneracy can be somewhat reduced but 
not eliminated by such a tuning, indicating the need for new terms in the action to reduce the nondegeneracy. 



1. INTRODUCTION 

Staggered fermions offer the possibility of do- 
ing unquenched calculations on current comput- 
ers with far less simulation time than Wilson type 
fermions. In their simplest form, however, they 
suffer from several well-known problems which 
must be addressed before they can be used effec- 
tively pj. One significant problem with ordinary 
staggered fermions is the large flavor nondegen- 
eracy, worst in the pion sector. At tree level, it 
arises from transitions between doubler quarks of 
different flavors induced by gluons of momentum 
ir dJl]. I n R e f- [B 1 Lepage showed how to turn 
the "fat-link" improvement of the MILC collabo- 
ration , which suppresses the coupling of quarks 
to these gluons, into a tree level 0(a 2 ) improved 
action by proper tuning of the coefficients and the 
inclusion of an additional term. 

Monte Carlo calculations showed a large reduc- 
tion in pion flavor nondegeneracy with this ac- 
tion Significant flavor breaking still remains, 
however, so further improvement is desirable for 
high precision calculations. In this work, we per- 
form a nonperturbative determination of gluonic 
corrections to the tadpole improved tree level a 2 
improved staggered action (called the "Asqtad" 
action by the MILC collaboration). We find a 
small improvement in pion flavor breaking, but 
not the large reduction that is still desirable. 
Therefore, incorporation of additional operators 
into the action is needed p| . 
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2. THE IMPROVED STAGGERED AC- 
TION 

We define a modified Asqtad (mAsqtad) action, 
built in two steps: First, from the naive fermionic 
action, one modifies the definition of the covariant 
derivative: 

D^i;(x) = V^(x + n)- V- M tl>(x - fj) (1) 

~lJf [(^) 3 ^ x + 3 ^ - ( U -^( x ~ 3M)] 
where is a fat link defined as 
V^ = g(l + co)L7 M 

£ u ±a u ±p u ±v u,u^ p u Ta 

-^E(^) 2C ^) 2 (2) 
o v 

(the indices in the sums are always different from 
H and among each other). The choice of the coef- 
ficients Ci — corresponds to the Asqtad action. 

Second, the fermion <p(x) is mapped into a 
scalar field x( x> ) by the relation 

4>(x) a a = £ {ltltltli l Y a x{x + A) (3) 
Ae[2 4 ] 

where a is the spin index, a is the flavor index, [2 4 ] 
is a 4D hypercube. Note that eq. (|J) is invertible 
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Figure 1. Three dimensional representation of 
the five dimensional space of coefficients and the 
points we considered. 



every value of j and k ^ jjh — 0.5). These points 
are represented in Fig. [y. This choice enables 
us to evaluate numerically the first (M| J and 
second derivative (M^) of the pion masses in 
respect to each coefficient of the action (keeping 
fixed the renormalization condition, i.e. the mass 
of the Goldstone pion). 

This work amounts to more than 10000 
fcrmionic inversions of the action and more than 
2000 fits of pion propagators; it was performed 
on the Fermilab QCD80 cluster. 



only if the fermion and the scalar live on different 
lattices, i.e. if the former lives on the blocked 
lattice of the latter . We use this prescription 
to build the 15-plet of pions of SU(4) flavor. 

Because of the remnant discrete flavor sym- 
metry there are only seven inequivalent pions. 
We identify the seven inequivalent pions by their 
SU(4) flavor structure 

C „,5 -,0_,5 ..3 .,5 ..1_,2 -.3.4 „,3 .,3 ( A \ 

?=7.77>7 7.7 7>1T>7.7 ( 4 ) 

Our goal is that of tuning the coefficients c 2 ; 
around zero to reduce the mass splitting among 
these pions. 

3. Computation 

Our computation was performed on 113 0(a 2 ) 
improved gauge configurations at /3 = 7.4 and 
Mo = 0.8629 on a 24 x 8 3 lattice. We preceded in 
the following way: 

We chose a finite set of points in the space of 
the coefficients. For each point, we performed 
a lattice computation of the Goldstone pion and 
we fine tuned the quark mass m in order to repro- 
duce a mass for the Goldstone pion, M 7 s , equal 
to (0.49 ± 0.01)a _1 (this is an arbitrary number). 
This fine tuning is required since the mass renor- 
malizes in different ways for the different choices 
of coefficients and we need to impose a physical 
renormalization condition. Then, for each point, 
we computed the whole spectrum of pions using 
the corresponding fine tuned value for the quark 
mass. 

We chose 91 points in the space of coefficients, 
namely Cj = 0, Cj = zszhdij, Cj = ihdij ± hSik for 



4. Results 

We expand the pion masses as functions of the 
coefficients in the action, in Taylor series up to 
second order 

AQ(ci) = Ml + J2 MiiPi + M liJ^ ( 5 ) 

i i,j 

and we define 

F 6 ( Ci ) = Mf(ci) - Afy (6) 

where M 7 s = const, because of our choice of 
renormalization condition. In Fig. |^ we report 
sections of the function F^. Each plot shows the 
value of F of all pions, £, when we vary one single 
coefficient. 

The first result that is already visible from the 
plots is that the spectrum is very mildly depen- 
dent on the action and there is no obvious direc- 
tion in the space of coefficients where the pion 
mass splitting gets reduced for all pions at once. 

We report here the values for the meson masses 
corresponding to the Asqtad action (for a fine 
tuned light quark mass of m — 0.0347a~ 1 ) 

M° 5 = 0.592 ± 0.003 ML 5 = 0.801 ± 0.017 

J° 7 7 

M°3 7 5 = 0.758 ± 0.006 M° lj2 = 1.021 ± 0.030 
M° 3 4 = 0.923 ±0.012 M° 3 = 1.092 ±0.032 
M°. = 0.977 ±0.015 

and their first derivatives 

M 7 o 7 s = (-0.022, -0.007, -0.008, 0.095, -0.036) 
M 7 3 7 s = (-0.040, -0.062, -0.020, 0.108, -0.007) 
ML 2 = (-0.005, -0.026, -0.016, 0.078, -0.035) 
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M^ 7 4 = (-0.036,-0.074,-0.029,0.103,-0.010) 
M' l3 = (+0.015, +0.006, -0.014, 0.071, -0.019) 
M^4 = (-0.022,-0.074,-0.034,0.067,-0.022) 

One should notice that the signs of the first 
derivatives are consistent for the different pions, 
but they are very small. They are of the same 
order of magnitude as the error on M^. This 
means that within a more than reasonable range 
(c, G [— 1,+1]) the pion's splitting does not vary 
more than 2a where a is its statistical error. 

We conclude that the flavor nondegeneracy of 
the tadpole improved tree level improved stag- 
gered action may be somewhat improved by a 
nonperturbative tuning of the coefficients, by per- 
haps 10% or less. However, for further dramatic 
reduction in flavor breaking, new terms must be 
added to the action, as in Ref. (!]]. 

This work was performed at Fermilab (U.S. De- 
partment of Energy Lab (operated by the Univer- 
sity Research Association, Inc.), under contract 
DE-AC02-76CHO3000. 
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Figure 2. F% as function of {cj}. 



